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Realistic methods to create vortices in spin-orbit-coupled Bose-Einstein condensates are discussed. 
It is shown that, contrary to common intuition, rotation of the trap containing a spin-orbit con- 
densate does not lead to an equilibrium state with static vortex structures, but gives rise instead 
to non-equilibrium behavior described by an intrinsically time-dependent Hamiltonian. We propose 
here the following alternative methods to induce thermodynamically stable static vortex configura- 
tions: (1) to rotate both the lasers and the anisotropic trap; and (2) to impose a synthetic Abelian 
field on top of synthetic spin-orbit interactions. Effective Hamiltonians for spin-orbit condensates 
under such perturbations are derived for most currently known realistic laser schemes that induce 
synthetic spin-orbit couplings. The Gross-Pit aevskii equation is solved for several experimentally 
relevant regimes. The new interesting effects include spatial separation of left- and right-moving 
spin-orbit condensates, the appearance of unusual vortex arrangements, and parity effects in vortex 
nucleation where the topological excitations are predicted to appear in pairs. All these phenomena 
are shown to be highly non-universal and depend strongly on a specific laser scheme and system 
parameters. 



I. INTRODUCTION 

Spin-orbit-coupled cold atoms represent a very new 
and quickly growing area of research that promises to 
host an even richer variety of exotic phenomena than 
solid-state spintronics [1 . Indeed, within just a few years 
of experimental research in the field, a number of exciting 
phenomena have already been observed [IHS] and there 
are clearly many more low-hanging fruits awaiting their 
experimental discovery. 

The key ideas underlying cold-atom spintronics - that 
studies particles with a synthetic spin degree of freedom 
coupled to their motion - grew out of the early theo- 
retical work by Juzeliunas et al. [6HT2], which showed 
that single-particle physics of atom-laser dressed states, 
where internal atomic states are coupled by position- 
dependent laser fields, can be described in terms of a 
non-Abelian vector potential acting on the dressed ex- 
citations. Later, it was demonstrated theoretically [13] 
that specific realizations of such laser configurations, in- 
cluding the early-proposed tripod scheme, give rise to 
spin-orbit-coupled Hamiltonians of Rashba-Dresselhaus 
type, familiar from solid-state semiconductor spintronics 
and that this "spintronics" description is a convenient al- 
ternative to the description in terms of the non-Abelian 
fields. Most importantly, it was quickly realized [14 that 
contrary to solid-state spintronics, where the underlying 
particles are bound to be electronic excitations, the syn- 
thetic spin- 1/2 degree of freedom in cold atoms can be 
carried by dressed spin-orbit-coupled bosons that were 
predicted to condense into a state, dubbed in Ref. 
a "spin-orbit coupled Bose-Einstein Condensate (BEC)." 
It was also shown [14] that multiple peaks in the time- 
of-flight expansion would be a smoking gun signature of 
such a new quantum state. Remarkably, this type of 
behavior was observed experimentally [2] by one of the 
authors shortly after. The specific laser setup used in 



Ref. ^ - that gives rise to an "Abelian" spin-orbit cou- 
pling (sometimes referred to as the "persistent-spin-helix 
symmetry point," [[T5l - [2T] where the Rashba and Dres- 
selhaus couplings are equal to each other) - was later 
analyzed in detail by Ho and collaborators in Ref. [22] . 
These experimental and theoretical successes have moti- 
vated other interesting theoretical proposals for realistic 
experimental schemes that can be used to create spin- 
orbit-coupled systems [23l [24] . Spin-orbit-coupled BECs 
have also been studied theoretically in Refs. [25H35] for 
different types of spin-orbit interactions and different in- 
ternal structures of bosons (pseudospin-1/2, spin-1 and 
spin-2 bosons). 

Among the obvious questions about the spin-orbit 
BECs is the physics of topological excitations - vortices 
- that play a central role in the physics of conventional 
BECs. This is subject of this paper, where we focus 
primarily on exploring experimentally-relevant methods 
that can be used to nucleate static vortex structures in 
spin-orbit BECs. In contrast to the conventional con- 
densates, the situation here is shown to be significantly 
more complicated as the vortex physics is obscured by 
the interplay of external perturbations intended to cre- 
ate them and the hyperfine structure underlying the syn- 
thetic spin-orbit-coupling setup. 

It is widely known and often taken for granted that 
rotating a Bose-Einstein condensate gives rise to the for- 
mation of vortices that arrange themselves into static 
vortex lattice structures. However, this picture is not in 
fact an obvious outcome of rotation, which represents a 
time- dependent perturbation due to a rotating anisotropic 
trap potential. The many subtleties involved in under- 
standing the fundamentals of the related phenomena are 
discussed in detail in the reviews by Leggett [36l[37], but 
the main conclusion is indeed that the physics of a one- 
component BEC confined to a spinning anisotropic trap 
can be mapped onto a statistical-mechanical problem of 
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the BEC with an effective time- independent Hamiltonian, 
i^eff = H — ujr ■ which describes the system in a ro- 
tating frame of reference (here, L is the orbital angular 
momentum operator and u)r is the frequency of rotation) . 

A naive expectation therefore is that to rotate an 
anisotropic trap would be a straightforward means to 
create vortex structures in spin-orbit-coupled BECs as 
well. However, this paper shows that this is generally 
not so and other, more sophisticated methods have to 
be involved in order to create static vortex structures. 
We show that the problem with rotation arises here 
because atoms are not influenced by the trapping po- 
tential only, but also by the lasers which create spin- 
orbit coupling in the first place. Therefore, if only the 
anisotropic potential rotates, it is in general impossible 
to choose a frame of reference where the Hamiltonian 
is time-independent, because the "spin-orbit coupling" 
lasers, stationary in the lab frame, are rotating in the ro- 
tating frame, generally resulting in non-trivial dynamics 
in any rotating frame. While there do exist rare degen- 
erate cases, where a unitary transformation that elimi- 
nates time-dependence from the non-interacting Hamil- 
tonian can be explicitly found, the interaction terms gen- 
erally become time-dependent under the unitary trans- 
formation, resulting again in a non-equilibrium problem. 
Hence, we argue that the residual time-dependence ap- 
pears to be an essential and unwelcome property of a 
spin-orbit-coupled BEC with a rotating anisotropic po- 
tential (at least for the realistic laser schemes currently 
known to us). We believe that while the specifics of time- 
evolution of rotating spin-orbit BECs are sensitive to de- 
tails of both the laser setup used and interactions, the 
typical scenario will involve non-universal dynamics that 
would inevitably lead to heating and destruction of the 
coherent state in contrast to the conventional BECs. 

It is therefore desirable to develop other 
experimentally-relevant methods to create vortices, 
like rotation or a magnetic field, for spin-orbit-coupled 
BECs. Two other ways suggested here and examined 
in detail are as follows: (i) to rotate both the lasers 
creating spin-orbit coupling and the trap, if the latter 
is anisotropic, or just the lasers for an isotropic trap 
(note that to rotate an isotropic trap has no meaning); 
(ii) To combine synthetic spin-orbit-couplings with a 
synthetic Abelian magnetic field. Theoretically, both 
methods are shown to give rise to interesting phenomena, 
including the appearance of sought-after static vortices 
and vortex lattices, parity effects in vortex nucleation, 
and real-space splitting of the spin-orbit BEC where the 
left- and right-moving parts are physically separated 
(an effect, which bears some similarity to the spin-Hall 
effect [38] |39j known in condensed matter spintronics) . 

Our paper is structured as follows: Sec. |TI| derives effec- 
tive Hamiltonians corresponding to a rotating trapping 
potential and/or rotating "spin-orbit lasers" for various 
spin-orbit-coupled laser schemes. In Sec. |III[ we solve the 
Gross-Pitaevskii equation to describe individual vortices 
and collective vortex structures for the laser scheme de- 



scribed in Ref. [7 with a rotating trap and Raman lasers. 
In SecjiVl we investigate vortex nucleation and other ef- 
fects associated with a synthetic magnetic field that can 
be imposed on top of the spin-orbit coupled system used 
in ^ by applying a spatially dependent Zeeman field. 



II. ROTATION IN SYSTEMS WITH 
ENGINEERED SPIN-ORBIT COUPLING 

In this section, we investigate the effect of rotation of 
an anisotropic trapping potential and/or spin-orbit lasers 
in three different laser schemes that have been proposed 
to create effective spin-orbit couplings. To distinguish 
between the different schemes, we will refer to the setup 
used in Ref. [4^ "M-scheme," the proposal described in 
Refs. [H [13] as "tripod-scheme," and the recent proposal 
of Ref. [2T as "4-level-scheme." 



A. M-scheme 

We first focus on the scheme used in recent experi- 
ment [4] and investigate the Hamiltonian for the case in 
which both trap and spin-orbit coupling lasers are ro- 
tating about the z-axis. The atoms in [4 are under the 
influence of three external sources: trapping potential, 
Raman lasers which create spin-orbit coupling and mag- 
netic field which creates Zeeman splitting (aligned along 
y-direction). If we wanted to get a time-independent 
Hamiltonian in the rotating frame we would have to ro- 
tate trapping potential, Raman lasers and magnetic field. 
To make things easier it is possible to change direction of 
the magnetic field to be along z-axis, which makes rota- 
tion of magnetic field about the z-axis unnecessary. If the 
change of the direction of magnetic field is accompanied 
by change in polarization of Raman lasers (the direction 
of lasers stays the same) the system is described by the 
same effective equations as in [4 . It is also important to 
note that, in the case of an isotropic trap, rotation of the 
trap does not have any effect and in that case rotating 
only the Raman lasers suffices. The stationary system is 
described by the following Hamiltonian (see methods in 
[4): 



Ha = 



2m 



V(v) 








(1) 



where k = — iV, V(r) is the trapping potential, 1 is the 
3x3 identity matrix, os^x^y^z are the 3x3 spin ma- 
trices, kL = V^tt/X, Q is the Raman coupling strength, 
ujz and ujq are the linear and quadratic Zeeman shifts, 
respectively. Here A is the wavelength and Aujl is the 
frequency difference of the two Raman beams used in 
the M-scheme. The Hamiltonian is written in the basis 
of hyperfine states {\mF = +1), = 0), \mF = — 1)} 
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which are quantized in z direction (direction of the ex- 
ternal magnetic field). 

When the trap and Raman lasers rotate with a con- 
stant frequency cj^ about the z-axis, the Hamiltonian 
^rot in the laboratory frame can be obtained from Eq. ([l]) 
using the following substitutions: 

V{x,y,z)^V{x(t),y(t),z) 



where 

x{t) = xcos{uJrt) + ysm{ujrt) 



(2) 



y{t) = ycos{uJrt) — xsm{ujrt) 
^3,x{t) = ^3,x cos{uJrt) + as^y sin(a;^t). 



(3) 



The Hamiltonian i^rot can be also written in a more com- 
pact form: 

where L is the orbital angular momentum operator and 
S is the spin operator, and and Sz are their z- 
components: Lz = h(xky — ykx)l^ Sz = has^z- 

The Hamiltonian Q is time-dependent in the labora- 
tory frame, but we show below that this time-dependence 
can be eliminated by a unitary transform. Recall that an 
arbitrary unitary transform, U{t)^ of the Hamiltonian H 
produces a new Hamiltonian, H\ as follows 



H' = UHU^ - ihU 



(5) 



We first go to the rotating frame of reference (rotating 
together with both the trap and the lasers) [40]: \iPrf) = 
U(t)\i)), where U{t) = exp[iujrt{Lz + Sz)/h]. Eq. ^ 
yields 



Hrf — Hq — uJr{Lz -\- Sz)j 



(6) 



where Hrf denotes the Hamiltonian in the rotating 
frame. The remaining time-dependence, arising from 
the oscillating Raman laser fields in Hq^ can be re- 
moved in the framework of the rotating wave approxi- 
mation. To obtain an effective description of the system 
in terms of two internal pseudo-spin states, we follow [4 
and choose the quadratic Zeeman shift huq to be large 
enough, so that the state \mz = 1) can be neglected. Us- 
ing the pseudo-spin- 1/2 labels for internal states, we get, 
1^) = \mz =0), 1^) = \mz = — 1). The final Hamiltonian 
can be expressed in the form used in Ref. [4] (a detailed 
derivation is much analogous to Ref. [4^ and is presented 
in appendix [A|) as follows. 



H 



RF,2 



2m 



V(v) 



Er 



kx^z + 77^x 

m 2 



hxIrkLyCTz 





— 5 



(7) 



where 1 is 2 x 2 unit matrix, cTx^y^z are 2x2 Pauli 
matrices and 5 = h{AujL — ^z) is a detuning from the 
Raman resonance. Since the resulting Hamiltonian is 
time-independent in the rotating frame, it leads to the 
appearance of stationary vortex structures studied below 
in Sec. [ml 

In the case where only the anisotropic trap is rotating, 
the Hamiltonian in the laboratory frame is given by ([T]), 
with V{x^y^z) V [x{t)^y{t)^ z) . Importantly, if we 
go to the rotating frame and make the rotating wave 
approximation (exactly as in the above), we are still left 
with a time-dependence (for details see appendix [B|) : 



^RF,2 — 



2m 



V(r) — UrLz 



(8) 



m 



'-kx{t)az 



where kx{t) = kx cos{uJrt) — ky sm{ujrt). 



B. Tripod scheme 

We now concentrate on the proposal described in 
Refs. [8 , 13 , which uses a so-called "tripod scheme," that 
consists of three degenerate ground states of an atom cou- 
pled to an excited state. The resulting energy spectrum 
includes two degenerate "dark" states and two "bright" 
states (one of the bright states is higher and the other is 
lower in energy with respect to degenerate dark states). 
In the strong coupling regime and within the adiabatic 
approximation, the energy difference between the dark 
and bright states is very large compared to other char- 
acteristic energies of the system. In this case, a cou- 
pling between the dark and bright states is negligible, 
and consequently if the atoms initially exist within the 
dark states subspace, they are expected to stay there for 
a long time. From now on, we use pseudo-spin- 1/2 nota- 
tions for the two degenerate dark states. 

The effective stationary Hamiltonian (projected onto 
the dark-state subspace) reads: 



H 



2m 



■ w{t) 



i - VoPx^y - ViPyCFz + ^QCFz, (9) 



where p = — ifiV, w{y) is a spin-independent part of the 
trapping potential (see appendix [C]for details), '^o and vi 
characterize the strength and type of spin-orbit coupling, 
and 5q is the effective Zeeman splitting. 1 is a 2 x 2 unit 
matrix, cTx.y.z are 2x2 Pauli matrices. 

We first investigate the case with both the trap and 
the spin-orbit lasers rotating. The derivation, presented 
in appendix [Cl leads to the following Hamiltonian in the 



4 



rotating frame 



exp[iu;^t(I/2/fi + (Jaj/2)], which gives: 



RF,2 



P 

2m 



• w{t) 



1 - VoPx^^^y ~ "^iPy^z 



SqCTz + mhWrivixcTz — Voy&y) 



(10) 



sm 

sm (/) cos (/) cos i 



sm cos (p cos f 
cos^ ^ cos^ 6 — sii 



where (j) = mvox/ cosO^ Sq 

( 



sin2^(6>/2) ) ]/^}/^' 6> is a constant. Let us note here 
that Ref. [41] previously considered the tripod scheme un- 
der rotation, but obtained shghtly different results (the 
spin angular momentum part {—uJrSz) was ignored in 
Ref. [H] ). 

Our result ([lO|, together with Eq. ^ for the M- 



scheme, clearly shows that the effect of rotation in sys- 
tems with synthetic spin-orbit interaction does not re- 
duce to just adding the —uj^Lz term for the Hamiltonian 
in the rotating frame, but also produces other position- 
dependent terms, which depend on a particular scheme. 

We now consider the tripod scheme with only the trap 
rotating. We first address the following question: if the 
trapping potential is time-dependent, can we get the ef- 
fective pseudo-spin Hamiltonian in the laboratory frame 
just by changing V V{t) in ([9|; or in other words, are 
we still allowed to restrict to the dark-state subspace if 
the external potential is time dependent? The answer is 
certainly "yes," if the trapping potential is the same for 
all three degenerate ground states (which is most often 
the case for optical trapping), because this kind of time- 
dependent potential does not couple the dark and bright 
states. 

In a general tripod scheme however, the trapping po- 
tential is not spin-independent (V^(r) = ^jVj{r)\j){j\^ 
Vi = V2 = w{t) and Vs = w{t) -\- S). To better under- 
stand this case, let us choose states {|1), |2), |3)} to be 
eigenstates of Sz (z-component of the total spin opera- 
tor). Then, the rotation of the trapping potential about 
the z-axis is described by: Vi = V2 = w'{Y^t) and V3 = 
w'{v,t)^5, where w'{Y,t) = e-^^-^^-/^^(r)e^^-^^-/^. 
We can therefore separate V{t) into a stationary 
spin-dependent term and a time-dependent but spin- 
independent term: V{r,t) = S\3){3\ + ^'(r, t) (|1)(1| + 
1 2) (2 1 + |3)(3|). Therefore the time-dependent part of 
trapping potential is spin-independent and it will not 
couple dark and bright states. With this, the tripod sys- 
tem with a trap rotating about the z-axis is described 
by 



H 



2m 



"^-VQpxCTy-VipyCTz^^Q^Z' (11) 



We now make the following transformation. U{t) 



H' 



where 



^ +^(r) -UJrLz 
2m 



ViPy{t)cFz{t) -^5Qcrz{t) 



1 - VoPx{t)ay{t) 

huJr 



(12) 



Px{t) = Px COs{uJrt) - Py sin{uJrt) , 
Py{t) = Py COs{uJrt) -\- Px sin(co'rt), 

d-y{t) = Gy cos{uJrt) — CTz sin(co'^t), 
d-z{t) = CFz cos{(jOrt) + Gy sin(cj^t). 



(13) 



The Hamiltonian (12) is generally time-dependent. How- 
ever in the case of Rashba coupling {vq = vi = v) and 
(5o = 0, this non-interacting part of the Hamiltonian be- 
comes static and reads. 



2m 



w{t) — UJrLz 
v{Px^y ^PyCTz) - 



1 



(14) 



C. 4- level scheme 

Here we study the 4-level-scheme [24 with a rotating 
trap. The stationary effective Hamiltonian (projected to 
the lowest energy states) is given by [24 : 



H 



21.2 



2m 



V{r) 



1 + a{axky - CTykx) 
A,, 



(15) 



P{cTxky + CFykx) + -^^z, 



where a and {3 denote strengths of Rashba and Dres- 
selhaus couplings respectively (in this scheme, a is fixed 
and j3 can be tuned), and is an effective Zeeman field. 
Per the same arguments as in the tripod scheme, we are 
allowed to simply replace V V{t) in (15) (if an exter- 
nal potential is time-dependent; note also, that the trap- 
ping potential here is spin independent). The rotating 



trap potential reads: F (r, t) = e 



rtL, 



rtL^/h 



We now make the following transformation: U{t) 
ex.-p[iuJrt{Lz/h-\- az/2)], which gives: 



V(r) - UJrLz 



1 + a{axky - CTykx) 



2m 

+ [a^y co^{2uJrt) + CTx sm(2ujrt)\ kx 
+ \&x cos{2ujrt) — cFy sin(2cj^t)] ky 



(16) 
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Again, this non-interacting part of the Hamiltonian is in 
general time-dependent, however for pure Rashba cou- 
pling, it becomes time-independent. 

Note that to get the full Hamiltonian in the rotat- 
ing frame, we must also include interactions between the 
bosons and apply to them the same transformations as 
in the non-interacting part above. If both the trap and 
spin-orbit lasers rotate, the corresponding unitary oper- 
ator, U{t) = ex.-p[iujrt{Lz + Sz)/h]^ describes a spatial 
rotation about the z-axis. If the bare interactions are 
rotationally-invariant, the interaction part of the Hamil- 
tonian does not change in the rotating frame. In con- 
trast to this result however, if only the trap is rotating, 
the interactions will generally acquire time-dependence 
as well (we have found a few very special cases - with se- 
rious constraints on the parameters of the system - where 
a unitary transform can be found that makes both the 
pure Rashba non-interacting part and interactions time- 
independent, but whether these degenerate cases can be 
realized experimentally remains unclear at this stage). 



III. CREATING VORTICES BY ROTATION 

In the previous section, we have shown that the Hamil- 
tonian for the M-scheme in the presence of a rotating trap 
and Raman lasers becomes time-independent in the ro- 
tating frame. In analogy with the physics of "ordinary" 
BEG under rotation, there will be thermal equilibration 
in the system and vortices will form in the condensate. 

Let us assume that huz ^ /J. {/J. is the chemical po- 
tential), which gives an effective 2D system, where the 
motion in z-direction is effectively frozen (this can be 
achieved by applying a ID optical lattice in z direction). 
We also assume the interaction part of the Hamiltonian 
to have the form: 



int 



d^r i ^Gip^ + ^G2pI + Gi2hPi ) > (17) 



where Gi, G2 and G12 are effective 2D interaction 
strengths and are related to 3D interaction strengths: 
Gi = Gf'^/iV2^k), G2 = Gfl{^U) and G12 = 
G?^/(\/2^l,), where I, = ^{ mUz). and pi are 
density operators for \i) states (normal ordering of 
the corresponding creation/annihilation operators is im- 
plied). 

We are interested in finding the ground state con- 
figuration of bosons in a rotating system described by 
( [7[l7| ). First, we have to make an assumption about the 
ground state and we assume below that (at the mean-field 
level) all atoms occupy the same single-particle state de- 
scribed by the spinor wave- function, (?/;^(r), '^^(r)) (we 
also call it condensate wave- function). The condensate 
wave-function satisfies the Gross-Pitaevskii (GP) equa- 



tions below: 



- ^ + ^ r - [Lz - hkLv) 

2m m ox ^ ^ 



+ NGM^+NGi2W 



2m ' ' " m dx 
5 + NG2\^Px? + NGi2\'>Pt? 



V{v) -ujr [Lz+hkiy- h) 



i>^ + ^i>^ (18) 



where N is the total number of particles and p is 
the Lagrange multiplier associated with the constraint 
J (i^r(|?/^^P + IV^iP) = 1 (it can be shown that p has a 
physical meaning of chemical potential [36 ). We solve 
the GP equations by using norm-preserving imaginary 
time propagation method (see for example Ref. [401 142]). 

We consider a trapping potential of the following form: 
V = ^muj'^{x^ + 7^^^)^ where uj and ^uj are trapping 
frequencies in the x and y direction. It is convenient to 
measure lengths in the units of the harmonic oscillator 
length, ao = \/h/ {muo) and energy in terms of hw. We 
introduce dimensionless position variable = r/ao- The 
corresponding "dimensionless GP equations" reads 



P^^ 



■ oj',{L', + kLv' - 1) + g2\i^i? + 5i2|^tl' 



^4- + IT^t. 



(19) 



where k'^ = ki^a^^ Vt' = Vt/{hjo)^ 5' = 5/{hw)^ uo'^ = 
cj^/cj, Lz = -i{x'dy' - y'd^j,'), gi = NGi/{hwal), g2 = 
NG2/{hujal) and gi2 = NG12/ ifu^jal). 

In simulations for the rotating system we consider ^^Rb 
atoms and we use the experimentally-relevant parame- 
ters: A = 804.1 nm, uo = 27r x 50 Hz and 7 = 1. These 
parameters give ao = ^JhJmuJ = 1.52 /im, k'j^ = 8.42. 
From now on we express length in units of ao (coordinates 
(x, y) in figures are also given in the units of ao). We per- 
formed simulations specifically for the rotation frequency 
uJr = 0.7 uj and for three different coupling strengths: no 
coupling ((^ = 0), weak coupling {Q = 2 E^)^ and strong 
coupling (yt = 10 El) {El = 35.4 hw). In simulations 
we choose gi = 1000, ^2 = 995, gi2 = 995. The ratio 
between ^1, ^2 and gi2 corresponds to interaction coef- 
ficients in ^^Rb (the interaction coefficients for ^^Rb in 
states {\F = l,m = 0),|F = l,m = —1)} is given in 
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Ref. [4]). In absence of rotation and for 11 = 0, (5 = our 
choice of ^i, g2 and gi2 produces cloud radius of 8.4 jiia. 
We also set 5 — hwr — 0- 

Without rotation and spin-orbit coupling, \X) and ||) 
components are miscible for our choice of interaction pa- 
rameters. In the case of rotation and no spin-orbit cou- 
pling there are several different phases depending on ujr 
and ratio of interaction coefficients [43 : triangular lat- 
tice, square lattice, stripe or double-core vortex lattice 
and vortex sheet. Since our Hamiltonian is almost equiv- 
alent to the Hamiltonian in Ref. [43] for (1 = and 
b — hujr = (there is a very small difference in interac- 
tion coefficients; the equivalence of non-interaction part 
of two systems is clear from (A3)) we reproduced results 
of Ref. [43J. 

The results for Q = are shown in Fig.[TJa), which dis- 
play the densities of the |t) and |^) components forming 
spatially-separated density stripes with lines of vortices 
along the minima of the density. As expected, our re- 
sults reproduce stripe vortex lattice phase described in 
Ref. |43l Note that for (1 = 0, the Hamiltonian ^ con- 
serves number of the |t) and |^) particles separately. We 
have chosen TV; {Ni J £r\ilji\'^). 

A weak spin-orbit coupling {Q = 2 El) (Fig.jljb)) does 
not appear to lead to any significant qualitative changes 
in the observed behavior: the densities of the |t) and 
\D components are still spatially separated and there are 
lines of vortices along the density minima of each com- 
ponent. 

A significant change comes in the strong-coupling 
regime: see the (1 = 10 El data shown in Fig. [ijc). 
The vortices arrange themselves in a lattice in |t) and 
\D components and densities of both components are al- 
most identical. This behavior can be understood from 
the following part of the Hamiltonian 



2m 



kL kx 



m 



(20) 



The spectrum of (20) for different ll's is shown in Fig. 
[2|a). For large it consists of two bands with an en- 
ergy separation much larger than all other characteristic 
energies of the system. Therefore, our system is "con- 
fined" to the lower band with a single minimum, which 
effectively makes it a single-component system. This ex- 
plains almost identical densities of the two components 
in Fig.jljc). 



IV. CREATING VORTICES BY 
SPATIALLY-DEPENDENT DETUNING 

A. The model 

Vortices in spin-orbit systems like [4 can be created 
without any rotation, but by imposing an additonal syn- 
thetic magnetic field. In [44 , it has been shown that a 
spatially-dependent detuning, ^, in the M-scheme results 
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FIG. 1: (color online) The density profiles for the rotat- 
ing spin-orbit-coupled BEG are shown. The first, second and 
third columns show density of |t) component den- 
sity of \D component (IV^^I^) and the total density {pr = 
I'^^l^ + respectively. Figure (a) shows results for Q = 

which are characterized by density stripes and lines of vortices 
in both components. The results for Q — 2 El (b) are qual- 
itatively similar to the Q = case. Figure (c) shows results 
for Q — 10 El] a vortex lattice is formed in both components 
and densities of the two components are almost identical. 



in a synthetic magnetic field, which creates vortices in 
the strong Raman coupling (Q) regime. Our goal is to 
investigate the same system for a wide range of Q (from 
weak to strong Raman coupling) and to see what kind of 
vortex structures it yields. 

The setup is described by the following effective Hamil- 
tonian (see j4ll44]): 



H 



2m 



V \1 



fi^kLkx 



m 



(21) 



We again assume strong confinement in the z direc- 
tion and describe interactions by equation (17). We are 



looking for the ground state in the same way as in the 
rotating case and following the same steps we get the 
"dimensionless GP equations": 



1 



/2 



5'iy') 



7 



2m 

<5'(y') 



7 + ^(^"+7Y^) 



^i + ^V'T- (22) 



Parameters ^7', (5^, /c^, ^i, ^12 are defined in tiie same 



way as in (19). 



B. Qualitative discussion 

To get a better understanding of the model, we inves- 
tigate Hamiltonian (21) in more detail. It is instructive 
to first focus on the following part of (21): 



- , h'^ kx ^ kL kx ^ 
li = — : i H (Jz 



2m 



rCTx 



m 



(23) 



We first assume that 5 is constant in space. In that case 
Hamiltonian (|23| can be easily diagonalized in the mo- 
mentum basis: W{kx)H\kx)U{kx) = 



E^{k^) 
E-{k^) 

The resulting spectrum consists of an upper(+) and 
lower(-) band, as shown in Fig. [2] The gap separat- 
ing the bands is large compared to other characteristic 
energies of the system and it is safe to assume that the 
condensate occupies only the states in the lower band. 
In Fig. |2|a), spectra for different coupling strengths Vt 
and ^ = are shown. For Vt < AEl^ the spectrum has 
two minima and BEC will involve states near both left 
and right minima. At ^7 = 4£^j^, there is a transition 
from a spectrum with two minima to a spectrum with 
one minimum, which changes the structure of the con- 
densate wave- function. I.e., for > AE^^ the BEC is 
expected to occupy only states with momentum around 
kx=^- 

The effect of detuning 5 in the low-1] regime is shown in 
Fig. ^h). We see that 8 shifts the energies and positions 
of the left and right minima. In the case of constant 5^ 
the BEC would occupy only the states around the global 
minimum (for example, the right minimum in Figj2|b)). 
Those cases have been tested experimentally in [4]7^ 

Now, consider a spatially-dependent 5{y). We will con- 
sider it to be a linear function of y: S{y) = So-\-Py^ which 
is the simplest and the most experimentally relevant 
regime. The interesting physics is evident from the fol- 
lowing arguments: for constant detuning, the spectrum 
around a minimum can be simply described by (we use 
dimensionless variables, see (19)) {kx — /^min)^/(2meff) + 

^rnin, whcrC 

^eff? ^min? -^min the cffcctivC maSS, 

position of the minimum, and the energy at the min- 
imum, respectively. Note that all these quantities de- 
pend on S. If ^ is ^/-dependent, the values of meff, 
^min, ^min will also bccomc spatially-dcpendent . Hence, 
the spectrum around the minimum can now be written 

as: [kx - kmin{y)] /(2meff) + E^in{y) which describes 
particles moving in an effective gauge field (A, <l>) = 
(^min(y), 0, 0, ^min(y)) with a Spatially- Varying effective 





(c) 



FIG. 2: (color online) The energy spectrum of H\ In (a) 
spectra for different Q (from Q = to Q = 6 El) and ^ = 
are shown (spectrum for Q = is at the top while spectrum 
for Q = 6 El is at the bottom). The effect of S in small Q 
regime is shown in (b) = 1 El, S = 0.5 El (solid blue line), 
^ = 1 El (dashed red line) and S = 2 El (dotted black line)). 
The effect of 6 in large Q regime is shown in (c) (Q = 16 El, 
6 = El (solid blue line), ^ = 1 El (dashed red line) and 
6 = 2 El (dotted black line)). 6 changes position and energy 
of the minimum. 



mass mQf[{y) [44]. The spatially-dependent vector poten- 
tial A induces an effective magnetic field (Bgff = V x A), 
which may lead to creation of vortices if strong enough. 
This approximation provides a good description of the 
system only if the particles at some point y have the 
momentum kx near the minimum. Our numerical simu- 
lations presented below indicate that this approximation 
in fact gives a very good qualitative description in a wide 
parameter range. 



We calculate parameters meff(^), /^^min(^), and ^min(^) 
by diagonalizing (23) for different y^s since 6 = S{y). The 
procedure of deriving effective equations for lower band 
for Hamiltonian (21) in high (single minimum) regime 
is described in Ref. [44]. Let us note however, that the 
method we use below to find the ground state is exact 
(in particular, we do not limit our system to lower band 
and we do not simplify interaction terms). 
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Results 



In simulations for a system with a spatially-dependent 
detuning S we use the same experimental parameters 
as in the simulations of a rotating system, which gives 
ao = ^yh/{muJ) — 1.52 /im and k'j^ = 8.42. We choose 
interaction parameters to be gi = 1600, g2 = 1593, 
9i2 = 1593 and constant part of detuning Sq = 0. 

The results for Q = 0, P = 4 hw/ao and 7 = 1 are 
shown in Fig. |3] and are straightforward to understand. 
In this case, we may write the Hamiltonian (21) as 



H, 



where = 

£ (k^ - 2kLk) + and V^{r) = V{r) + S{y)/2, 

Vi{r) = V(r) - S{y)/2. We see that motion of |t) and |;) 
particles is decoupled in H and that they experience dif- 
ferent potentials V^{t), V;(r). Detuning gradient /3 shifts 
the minima of V^{r) (V;(r)) for yo = /3/(2mco'^7^) in the 
positive/negative ^-direction and therefore, the centers 
of the It) and ||) densities are shifted from the origin by 
±?/o Cthe origin is located in the minimum of ^(r)), see 
Fig. [i^b). Also, it is clear from H and Fig. [ija) that the 
momentum distribution of |t) particles will be cen- 
tered around k = (— /cl,0) (k = (/cl,0)), see Fig. [sj^c). 
The effect of repulsive interactions between the particles 
with different spins is clearly seen (the overlap between 
It) and It) densities is quite small). 
If we introduce a finite Q, the Hamiltonian becomes: 



H 



Q/2 I ' "^^^ corresponding (^-term creates cou- 
pling between |t) and |t) particles. If S = const = and 
ft is small, the states around the left (right) minimum in 
the spectrum in Fig. [2|a) still consist mainly of the |t) 
(It)) particles, but there is also some admixture of the 
component with the opposite spin, which grows with Q. 
It means that '^^(r) (?/^|(r)) will mainly consist of states 
with momentum around the left (right) minimum, but 
also of states around the right (left) minimum. We can 
therefore write: 



(24) 



where and V^|L(r) consist only of states with mo- 

menta around left peak, while ip^R{r) and ipiR{r) consist 
only of the states with momenta around right peak of mo- 
mentum distribution. We therefore call (?/;^L(r), '^|L(r)) 
and (^/^^^^(r), ?/;|i^(r)) left and right wave- function. In the 
spatially-dependent detuning case it may happen that 
momentum distribution is separated in two peaks (i.e., 
there exist "left"- and "right-movers") even for Q > 4 El 
(see for example Fig. [tJ^c)). In that case also the notion 
of left and right wave-function applies. 

To investigate the effect of which couples |t) and 
It) states, we consider the regime with Q = 3 El and 
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FIG. 3: (color online) The figure shows results for r2 = 0, 
P = 4 hu/ao and 7 = 1. In (a) the total density is shown. 
The shape of the density is determined by spatially-dependent 
detuning, which shifts the densities of |t) and |t) particles (b). 
Momentum distribution of |t) and |t) components is shown 
in (c). 



P = 8 huj/ao (Fig. [4|. The total density Pt{'^) is shown 
in Fig. 4]^ a) and there is a characteristic series of minima 
along the x-direction at 7/ = 0, which come from vortices 
in the ?/^^ and wave- functions, see Fig.|4|b), which are 
positioned along x and near y = 0. We have checked that 
the phase winding around zero density points of |?/^|p and 
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IV^^P is — 27r. Since vortices in ||) and |t) components 
are slightly displaced from ^ = 0, the density at minima 
in pT are close to, but not exactly equal to zero. To 
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kt(k)r 




the amplitude of ip^R is considerably smaller than the 
amplitude of ifj^L'- J c^^^l'^t^P = 0.05 and / (i^r|V^^Lp = 
0.45). The momentum distribution in Fig. ^c) shows 
that the wave-packet, ^|J^L^ has an average momentum of 
^left = —0.8 kL and ip^R has an average momentum of 
bright = 0.8 Since ip^ is a superposition of the left- 
and right-movers, ip^ = V^^L + V^t^' appearance of the 
line of vortices at overlapping region is expected. The 
separation of vortices d is then simply given by (/Cright — 
^ieft)<^ = 27r or (i = (27r)/(A:right — ^left)- The analytical 
expression for d fits perfectly well to our numerical data. 
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FIG. 4: (color online) The figure shows results for Q = 3 El, 
(3 — 8 hw/ao and 7 = 1. In (a) the total density is shown. 
The series of minima at = comes from vortices in |t) and 
\\) wavef unctions (b). Momentum distribution of |t) and |^) 
components is shown in (c). 

explain the existence of the line of vortices in the W) 
and 1^) components, we examine the left and right wave- 
functions. Fig. p[a) displays IV^^lP and I'^^i^P (note that 



(b) 

FIG. 5: (color online) The figure (a) shows iV^t-f^l^ IV^T-rI^ 
the relative amplitude of which is given by / d^r\il:^L\^ — 
0.45 and / d^r\ij;^R\^ = 0.05 for the parameters r2 = 3 El, 
/3 = 8 hjj/aQ and 7=1. The superposition of ij^^L and iI^^r, 
— ^^j^ -\- ip^R^ produces vortices in ip^. The density of 
left- and right-moving particles {pL — \'P-\l\^ + IV^^lI^, Pr — 
iV^t-^l^ + iV^^i?!^) particles is shown in (b). 

To explain the density profile and momentum dis- 
tribution, it is useful to consider an effective gauge- 
field picture. The effective gauge field, (A, ^) = 
(^min(y), 0, 0, £^min(y)), cau be Calculated by diagonal- 
izing H' . As discussed earlier, we may approximate 
the low-energy band physics by the following Hamilto- 
nian (we use again the dimensionless variables, where the 
lengths are measured in terms of ao and the wave- vectors. 



10 



k, in terms of 1/ao): 



Hegf 



1 



2TOeff(y) 



[k.-A{y)y 



My) + V{r), (25) 



-tV). 



For Q > AEl there is a sin- 



where V{y) — \{x^ 
gle local minimum in lower band of the Hamiltonian ( 23 ) 
spectrum for any 5. For Vt < AE^ the spectrum has two 
minima for (5 = 0, however when 5 becomes large enough 
the spectrum has a single local minimum (Fig. |2|b)). 
The spectrum around each local minimum can be ap- 
proximated by the form given in ( [25| ) , and therefore there 
will be AL{y)^ ^L{y), ^eff,L(^) corresponding to the left 
minimum and Aii{y)^ ^ii{y), ^eff,R(^) corresponding to 
the right minimum of the spectrum. Left-movers feel the 
"left gauge field" (^/^(y), 0, 0, <I>l(^)) while right-movers 
feel the "right gauge field" {AR{y),0,0,^R{y)). 

To get the effective potential in y direction acting on 
left- and right-movers we define: leff,L(^) = ^Liv) + 



|7^y^ Vei^Mv) = ^niy) + h'^y'^' in Fig. ^ we 

show ^L/R{y)^ '^eff,L/R(^), ^L/i?(^) and l/meff,L/R(^) 

for 1] = 3 El and /3 = 8 hu/ao. V'eff,L/R have minima 
at yo^R/L — i3.2 which explains the total density profile 
(Fig.pfa)) which has maxima at = ±3.2. The position 
of two peaks in momentum distribution in Fig. |4jc) can 
be understood as follows: for particles positioned near 
the minimum of 14ff,L in Fig- [sj^b), it is energetically fa- 
vorable to have the f-component of momentum approxi- 
mately equal to A(?/o,l) and the ^-component near zero. 
Fig.[6|c) shows that A(?/o,l) ^ —0.79 /cl, while from Fig. 
|4](c), we see that the momentum distribution is centered 
around kx — —0.80 /c^. The same explanation applies for 
the momentum distribution of right-movers. 

To investigate the regime with a single minimum in 
the spectrum (Q. > AE^) we did calculations for parame- 
ters: O = bEL, P = l2huj/ao and 7 = 1 (Fig. [7|). In this 
"single-minimum" case one might expect momentum dis- 
tribution to be concentrated around a single point as was 
observed in Ref. [¥. However, in spatially-dependent- 
detuning case this is not necessarily true: the momen- 
tum distribution (Fig. [tJ^c)) shows two peaks around 
kx = ±0.55 kL. Also, the total density (Fig. [7|^a)) has a 
characteristic series of minima along ^ = line which 
come from vortices in the tjj^ and t/j^ wave-functions 
(Fig. [rj^b)) created in the overlapping region of left- and 
right-movers. The results can again be explained by the 
effective gauge field. The effective potential in ^-direction 
Kff(^) = ^{y) ± ^7^^^ (Fig. [sj^a)) has two minima at 
yo,R/L = ±3.4 which explains the density distribution 
which has maxima at = ±3.3. Also, equation (23) tells 



us it is energetically favourable for particles near the left 
(right) minimum of V^ff(^) to have momentum around 
A(yo,L) = -0.56 kL {A{yo,R) = 0.56 fc^) (Fig. gc)) 
which explains momentum distribution. We also note 
that in Fig. [sjc) A{y) has a large gradient and therefore 
magnetic field (^eff ^ dA/dy) is strong around y = 
which may serve as an alternative explanation of line of 
vortices appearing in Fig. [7|^a). 

We now study the system with strong Raman cou- 





E 0.5 



FIG. 6: (color online) The figure shows the scalar potential 
^{y) (a), the effective trapping potential in y-direction Vef[{y) 
(b), vector potential A{y) (c) and inverse of the effective mass 
(d) for Q = 3 El J 13 = 8 huo/ao and 7=1. Values correspond- 
ing to the left minimum of the spectrum are represented by a 
solid red line while the values corresponding to the left min- 
imum of the spectrum are represented by a dashed blue line 
(see the text for details). 



pling ft and weak detuning gradient f3 (i.e. /3 is not 
large enough to produce spatial separation of a cloud 
along y as in previous cases). Results for ^7 = 10 El, 
(3 = 12 hhj/ao are shown in Fig. [o] and can be explained 
by the associated effective gauge field shown in Fig. [To] 
The total density (Fig.jo^a)) and IV^^p, (Fig.jgjb)) 
show the existence of a vortex in the centre of the cloud. 
The vortex appears only for strong enough effective mag- 
netic field which is tuned by changing p. We define the 
effective magnetic field Beff = V x A(?/) and in our 
case, (A = (A(?/), 0, 0)), Beff = -^^z. The magnetic 
field points in the z direction, depends on y, and is con- 
stant along X. We also note that since mQf[{y) ^ 1 (Fig. 
[Toj^d)), the effective equations will differ from those for 
an ordinary charged particle in a magnetic field 5eff(^)^- 
The vector potential A{y) and the effective magnetic field 
Beff(y) are shown in Fig. 10'b,c). 



It is useful to know the critical field needed for vor- 
tex creation and we may get a crude estimate by us- 
ing the equation for critical magnetic field of a single- 
component 2D gas in the Thomas-Fermi limit: Be = 
4(ao/i^)^ln (0.888(i?/ao)^), where R is the Thomas- 
Fermi radius of the cloud |45]. We take R = 6.5 ao (the 
size of our cloud), which gives ~ 0.35. It is important 
to notice that larger number of particles or stronger in- 
teractions increase R, which lowers the critical field {Be 
decreases with increasing R). To find Be, we did simu- 
lations for 1] = 10 El, 7 = 1 and for different values of 
P (which controls the strength of the effective magnetic 
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FIG. 8: (color online) The figure shows the scalar potential 
^(y) (a), the effective trapping potential in y-direction Vef[{y) 
(b), vector potential A{y) (c) and inverse of the effective mass 
(d) for Q = 5 El, (3 = 12 huj/ao and 7=1. 
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FIG. 7: (color online) The figure shows results for Q = 5 El, 
(3 = 12 hu/ao and 7=1. In (a) the total density is shown. 
The series of minima at 2/ = comes from vortices in |t) and 
\D wavefunctions (b). Momentum distribution of |t) and |^) 
components is shown in (c). 



field). We found that the vortices start to appear for a 
critical effective magnetic field Be ~ 0.34, which is very 
close to our estimate presented above. 

If the effective field is strong enough, a vortex "lat- 
tice" is formed, as shown in Fig. [TTJ which corresponds 
to = 10 El, /3 = 40 huj/ao and 7 = 1.85. From the 
figure, we see that vortices are concentrated along the 



X-axis and around ^ = 0. This is because 5eff(y) is not 
homogeneous, i.e. the field is strongest at y = and 
it weakens with increasing \y\. We had to increase trap- 
ping strength in the y direction (7 = 1.85) because scalar 
potential ^{y) separates the clouds (e.g. see Fig. [sj^a)) 
and for a weaker trapping strength, the effective poten- 
tial would have two minima (it would look like effective 
potential in Fig. |8|b)). 

The most interesting regime is the one in which left 
and right moving phases (?/;^L(r), '^t^(r)) and [iIj^r{t), 
i^tRi^)) spatially separated along y direction and 
there is a vortex (or vortices) in each phase in addition to 
a vortex line. This requires double minimum structure of 
the effective potential in y direction T4ff(^), which sep- 
arates the phases and strong enough effective magnetic 
field in each phase to create additional vortices, which 
tend to appear in pairs (i.e. the number of vortices is 
equal in both phases which is a consequence of the fact 
that in our simulations the effective gauge field is sym- 
metric with respect to reflection about ^ = line and 
interactions are almost spin-independent). 

In Fig. [l2|a), we show results for a; = 27r x 10 Hz, Q = 
4 El^ P = 20 hw/ao (ao = ^/h/ {muj)). By choosing uj = 
27r X 10 Hz, parameter /c^ in dimensionless GP equations 
(22) becomes /c^ = 18.83, while interaction coefficients 
stay the same (^1 = 1600, ^2 = 1593 and gi2 = 1593). 
Having larger /c^ means we can create stronger effective 
magnetic field. We increased the trapping frequency in 
the ^/-direction (7 = 1.3) to bring two phases closer to 
y = 0, where the effective field is stronger (to counter 
the effective scalar potential ^{y), which separates the 
phases). In Fig. 12 'b) we show results for a; = 27rx 10 Hz, 
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FIG. 9: (color online) The figure shows results for Q = 10 El, 
/3 = 12 h(jj/ao and 7=1. In (a) the total density is shown, 
while (b) and (c) show densities of |t) and \D (c) components. 
The vortex in the center appears for strong enough effective 
magnetic field. 



n = 10 El, P = 150 huj/ao and 7 = 2.75. Here the left 
and right phases are completely separated in space and 
the effective magnetic field is strong enough to produce 
multiple vortices in each phase. Also, it is clear that the 
vortices are not located in centers of two phases, but are 
positioned closer to y = which is expected because the 
field is stronger near y = 0. 

It is important to discuss the means of experimen- 
tally observing results we presented. We concentrate 
on the time-of-flight imaging, which is widely used to 
probe cold- atoms systems. The time-of-flight picture 
here will be determined by the underlying momentum 
distribution of particles. If this momentum distribu- 
tion consists of two separated peaks, the initial cloud 
will strongly separate during expansion (see for exam- 
ple [SI |4j). We note that due to the transformation 
'0|(r) = '^|(r)e~*^^^, = '0^(r)e*^^^ used when de- 

riving Hamiltonian ( [2T| , the real momentum distribu- 
tion of It) particles will in fact be shifted by kL with 
respect to the momentum distribution shown in figures 
and the momentum distribution of ||) particles is shifted 
by —kL (see [4 ). In the case of 1] = 0, both |t) and 
1^) particles will have zero average momentum, which 
means both components of the condensate will expand, 
while the position of centre of mass will be stationary 
during time-of-flight. For Q = 3 El and /3 = 8 hu/ao, 
we expect four separated clouds to be seen in the time-of- 
fiight: since the real momentum distributions of |t) and 
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FIG. 10: (color online) The figure shows the effective trapping 
potential in y-direction Vefi{y) (a), vector potential A{y) (b), 
the effective magnetic field Beff (c) and inverse of the effective 
mass (d) for Q = 10 El, /3 = 12 huj/ao and 7=1. 



5l 

>> ol 

5I 





X 



FIG. 11: (color online) The figure shows the total density for 
Q = 10 El, 13 = 40 huj/ao and 7 = 1.85. 



\D particles are shifted by kL and — /cl, there will be two 
clouds of It) particles with average momenta of 0.2 kL 
(larger cloud) and 1.8 /cl (smaller cloud) and two clouds 
of \D particles with average momenta of —0.2 kL (larger 
cloud) and —1.8 kL (smaller cloud). It is important to 
notice that the vortex line will not be easily visible in 
those images, because it exists only due to the the over- 
lap of the wave-packets with different average momenta. 
During the time-of-flight, two wave-packets {ip^,L,i^t,R) 
or {i^i,L,i^i,R) separate, which means that they do not 
overlap any more and there is no clear vortex line present. 



For the case in Fig. 12, the vortices in each phase will be 
visible since they are not a result of overlapping the left- 
and right-moving condensates. 
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FIG. 12: (color online) Figures show separated left and right 
phases with vortices in each phase. Trapping frequency is 
cj = 27r X 10 Hz. Figure (a) shows total density for Q = 4 El, 
/3 = 20 huu/ao and 7 = 1.3. Figure (b) shows total density for 
Q = 10 El, 13 = 150 hw/ao and 7 2.75. 



V. CONCLUSION 

In this paper, we have investigated realistic experi- 
mental methods that can be used to create vortex exci- 
tations in spin-orbit-coupled Bose-Einstein condensates. 
The main conclusion of the work is that due to a com- 
plicated interplay between effects associated with the ap- 
plied laser fields and rotation, the resulting state of the 
spin-orbit BEG under additional perturbations is highly 
non-universal and depends strongly on the system pa- 
rameters and specific laser schemes. In particular, we 
argued that a spin-orbit BEG under rotation of the trap 
alone does not achieve a thermodynamically stable state 
at all, but acquires a complicated non-equilibrium dy- 
namics that eventually leads to heating and the destruc- 
tion of the condensate. 

We have also suggested two alternative experimen- 
tal methods to mimic an Abelian "orbital" magnetic 
field that involve either rotation of the entire experimen- 
tal setup, or a spatially-dependent detuning. We per- 
formed numerical simulations of the resulting thermody- 
namically stable density distributions, focusing mostly on 
the M-scheme that has already been realized experimen- 
tally. This scheme gives rise to an "Abelian" spin-orbit- 
coupling with a well-understood ground state that we 
used as a basis of our numerical simulations that showed 
topological excitations above the ground state. We ex- 
pect that the predicted vortex configurations, in particu- 
lar vortices appearing in pairs in the spatially-separated 
left- and right-moving regions, would be straightforward 
to observe experimentally, as all necessary ingredients are 
already experimentally available. 

Finally, we would like to mention that to this point 
only an "Abelian" spin-orbit-coupling scheme has been 
actually realized in experiment [4 and we mostly focused 
here on vortex topological excitations in such systems 



with a well-understood ground state. What remains of 
great interest of course is an expermiental realization of 
a truly "non- Abelian" spin-orbit interaction (either of 
pure Rashba or Dresselhaus type or a non-equal mix- 
ture of those) , which can be achieved using laser schemes 
described in Sees. |ll^, [Hp, and Refs. [131 M- Note 
that it was argued theoretically [14 that in the Rashba- 
Dresselhaus system with single-particle dispersion of the 
double-well type, a fragmented condensed state [46l [47] 
can be selected by energetics for repulsive interactions 
that do not break the underlynig Kramers symmetry. 
This state arises because repulsive interactions in the real 
space tend to localize particles in the dual momentum 
space per the fundamental Heisenberg uncertainty prin- 
ciple. This robust argument together with the protection 
provided by Kramers symmetry and momentum conser- 
vation (modulo finite-size effects due to the trap) suggest 
that the long-sought-after fragmented BEG (which takes 
the form of a many-body Schrodinger's cat-state in this 
case [14 ) is more stable in spin-orbit-coupled systems 
than that in BEGs confined to real-space double-well po- 
tentials and hence can be observed experimentally. Topo- 
logical excitations above this exotic ground state are ex- 
pected to also be of exotic nature and may potentially re- 
alize much of the exciting physics discussed in the context 
of multi-component superconductors [48l [49]. Finally, 
the nature of the ground state and topological excita- 
tions above it in the pure bosonic Rashba model remain 
of great interest as well. Depending on the interaction pa- 
rameters, this model with a continuous ring of minima on 
a circle in momentum space, may potentially host topo- 
logical BEGs, spontaneous symmetry-broken states [33], 
and exotic Bose-liquid states [50] , where strong quantum 
fluctuations prohibit order even at zero temperature. 
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Note added - After this work was completed, we be- 
came aware of two papers which study spin-orbit-coupled 
BEGs under rotation [511 [52] . The fundamental assump- 
tions in these papers are qualitatively different from our 
theory, in that Refs. [51, 52j start with an effective spin- 
orbit-coupled Hamiltonian and assume that it remains 
stationary under rotation. This is in contrast to our the- 
ory, where we consider realistic experimental schemes, 
where rotation is shown to lead to a different descrip- 
tion. 
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Appendix A: M-scheme with rotating trap and 
spin-orbit lasers 



The Hamiltonian in the rotating frame (IgI) is: 



Hrf = 



2m 



V(v)-UrL, i 



(Al) 





hu^ 



The Hamiltonian becomes time-independent if we trans- 
fer to the rotating-wave frame and if we do the rotating- 
wave approximation: 



Hrf 



2m 



V{t) — UOrLz 



^s^hujq 

-6j 



(A2) 

where 5 = h{AujL — ^z)- We set quadratic Zeeman shift 
hujq to be much greater than Q and S so we may restrict 
to the subspace spanned by {\mz = 0), Im^ = —1)}: 



H 



RF,2 



21.2 



2m 



V(r) — uJrLz 



--a,sm(2fe^x)+^Q 



(A3) 



where 1 is 2 x 2 unit matrix and cfx^y^z are 2x2 Pauli 
matrices. Since there are effecively two internal degrees 
of freedom we introduce pseudospin-1/2 notation, i.e. 
we define |t) = = 0), |;) = = -1). We fol- 
low the steps in [4 and make transformation: '^|(r) = 
V^^(r)e-^^^^, V^|(r) = V^;(r)e^^^^, where (V^^(r), V^;(r)) is 
a spinor wavefunction on which Hamiltonian (A3) acts. 
The Hamiltonian then becomes: 



F,2 



2m 



V{r) - ujrL^ + El 



(A4) 



m 



where El = h^k\/2m. We can drop E^^ term by simply 
renormalizing the energy. 



Appendix B: M-scheme with rotating trap 

The Hamiltonian H[.^^ describing M-scheme with ro- 
tating trap in the laboratory frame is: 



H' 



21^2 



V{x{t),y{t),z) - UJrLz 



2m 

'h{—UJz + UJq) 





hwz 



(Bl) 



-h \/2^lcT2>,x cos{2kLX + Auj^t)^ 



where x{t) is defined in (|3|. After transfering to the 
rotating frame {U{t) = ex.-p[iuJrt{Lz-\-Sz)/h]) and making 
the rotating wave approximation the Hamiltonian is: 



Hrf = 



2m 



V(r) — UJrLz 



1 - hUrOS z 



^35/2 ^hujq 
6/2 
-6/2j 



(B2) 



n 

T2 



G'^^x COs{2kLx' {t) + UJrt) 

as^y sm{2kLx\t) + cj^t), 



where x'{t) = xcos{uJrt) — ysm{ujrt). We may again 
neglect state \mz = 1) assuming ujq ». To get 
the Hamiltonian in a more familiar spin-orbit-coupling 
form we make the following transformation: '0|(r) = 

V^t(r)e-^^^^'^^\ V^|(r) = ^ij^(^r)e'^^'^^^^^+'^^^ , which gives: 



Hrf, 2 — 



2m 

h^kL 
m 



V{v) - UJrLz + El 



(B3) 



tCTx 



where kxit) = kxCos{uJrt) — ky sm{ujrt). We can drop 
ElI term by renormalizing the energy. 



Appendix C: Tripod scheme with rotating trap and 
spin-orbit lasers 

The original Hamiltonian for the tripod scheme (sta- 
tionary system) is (see [13 ): 

Ho = ^i + V{r) + Ha-i, (CI) 
2m 

where V{t) = V^(r)|j)(j| is spin dependent trap- 
ping potential, atom-laser interaction Ha-i = A|0)(0| — 
((^i|0)(l| + 1^2|0)(2| + ^^3|0)(3| -h H.c.) and i is 4 X 4 
unit matrix. A is detuning from resonance and ^1,2,3 
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are Rabi frequencies: fti{r) = fls'mO cos{mVax)e'^^'^^y ^ 
(^2(r) = I]sinl9sin(mva2;)e^^^^^, l]i(r) = ncosO, 
where 0^ Va and v^j are constants (see [13] for 
details). If we start rotating spin-orbit lasers in 
the laboratory, atom-laser interaction part of the 
Hamiltonian becomes e-^^-^(^-+'^-)/^^^_^e^^-^(^-+'^-)/^. 
If the trap rotates, trapping potential becomes 
^-iu;r^tiL,+s,)/hY^iujr^tiL,+s,)/h^ Therefore, we can write 
the Hamiltonian of the rotating system as: 



(C2) 



The Hamiltonian in the rotating frame is then: Hrf = 
Ho — uJr{Lz-\-Sz)' Since Hrf is time-independent we can 
use exactly the same procedure for getting the effective 
spin-orbit coupling described in |8l[T3^, i.e. we project 
the Hamiltonian to the dark states subspace. Here we 



assume that three degenerate hyperfine groundstates are 
part of F=l manifold (for example the ground state of 
^^Rb) and that they are eigenstates of Sz- This gives 
us the precise form of Sz operator. As in [13 we take 
Vi = V2 = w{r) and Vs = w{r) + S. After projecting to 
dark states we get: 



H 



RF,2 



2m 



-^w{r)-Lz 



i - VoPxCTy - ViPy&z 



f s\v? (j) sin (j) cos (j) cos 

^ \ sin 6 cos 6 cos cos^ cos^ 6 — sin^ / 



(C3) 



where 5o = sin^ e{5 - + (-^f]/2}/2. 
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